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SOME  APPROXIMATE  BOLTZMANN  EQU  TION  SOLUTION  METHODS 
FOR  THE  STUDY  OF  INTENSE  MICROWAVE  PULSE  PROPAGATION 


I.  INTRODUCTION 

The  propagation  through  the  atmosphere  o£  short  pulse  high  power 
microwave  (HPM)  radiations  is  of  interest  for  many  applications.  One  of  the 
important  aspects  of  the  propagation,  however,  is  the  attenuation  of  the  pulse 
as  it  traverses  the  atmosphere.  The  attenuation  results  from  the  absorption 
of  the  radiation  by  the  air  plasma  electrons,  which  are  generated  by  the  air 
breakdown  action  of  the  high  power  microwave.  The  microwave  radiation  is  also 
attenuated  through  absorption  by  the  air  molecules  and  by  scattering  processes 
with  molecules  and  particulates. 

To  calculate  the  attenuation  of  the  HPM  pulse  one  must  solve  Maxwell's 
equations  in  conjunction  with  the  radiative  transfer  equation.  The 
attenuation  coefficient,  due  to  the  absorption  by  the  air  plasma,  is  related 
to  the  plasma  conductivity  or  more  specifically  it  depends  on  the  electron 
density  and  the  collision  frequency  for  the  momentum  transfer.  The  electron 
density,  on  the  other  hand,  depends  on  the  avalanche  ionization  rate  which  in 
turn  depends  on  the  atomic  cross  sections  and  the  appropriate  electron 
velocity  distributions. 

The  relevant  quantities  for  the  HPM  attenuation  calculations,  however, 
can  be  obtained  if  an  ab  initio  approach  is  considered,  where  one  solves  the 
Boltzmann  equation  in  conjunction  with  Maxwell's  equations  and  the  radiative 
transfer  model. 

To  obtain  reasonable  estimates  for  the  microwave  attenuation,  one  may  use 
experimental  values  for  the  ionization  and  momentum  transfer  rates  and  w  . 
However,  in  the  region  of  interest,  i.e.  high  power  and  short  pulse,  there 
exist  no  experimental  data.  Therefore,  the  Boltzmann  equation  approach  is 
necessary.  The  formulation  of  a  Boltzmann  equation  requires  a  set  of 
appropriate  atomic  cross  sections  which  represent  the  state  of  the  art  and  are 
reliable.  The  appropriate  cross  sections  are  given  in  tabular  form  at  the  end 
of  this  report  based  on  a  recent  compilation  by  Ali^. 

Manuscript  approved  July  23.  1986. 


The  Boltzmann  equation  approach  in  obtaining  the  electron  velocity 

distribution  has  been  formulated  in  general  by  many  authors  and  its  solution 

2  3  A 

using  spherical  harmonics  is  veil  understood'1’  ’  .  Many  treatments  of  the 

Boltzmann  equation,  applied  to  breakdown  in  gases  and  swarm  experiments,  have 

utilized  the  steady  state  approach  where  ionization  effects  are  not 

considered.  This  approach  has  also  been  used  previously  as  means  to  obtain 

electron  atom  collision  cross  sections^. 

The  effect  of  ionization  for  high  E/N  (E  is  the  electric  field  and  N  the 

gas  density)  has  been  treated  for  swarm  experiments  by  Taniguchi,  et.  al.,6 
7  8 

Brunet  and  Vincent  and  Yoshida,  et.  al  .  The  electron  velocity  distribution 

o 

in  gas  breakdown  by  an  rf  pulse  has  been  treated  by  Morez  .  However, 
Boltzmann  equation  analyses  for  very  high  rf  fields  with  short  pulse  durations 
have  had  very  little  attention.  Recently,  Pitchford^* ^  has  considered  the 
time  dependent  Boltzmann  approach  in  the  breakdown  of  nitrogen. 

In  the  following  sections,  we  describe  some  theoretical  formulations 
which  provide  a  basis  for  generating  approximate  solutions  of  the  Boltzmann 
equation  suitable  for  the  study  of  short  intense  microwave  pulses  interacting 
with  a  background  atmosphere,  in  various  parameter  regimes.  It  is  not 
intended  that  all  of  these  formulations  be  developed  in  full  detail  here; 
rather  we  shall  simply  outline  calculation  methods  being  considered  in  order 
to  provide  a  platform  for  discussion  of  the  advantages  and  disadvantages  of 
the  various  approaches.  One  example,  appropriate  in  the  collisional  regime, 
is  developed  in  detail  and  implemented  on  a  sample  problem. 

II.  TIMESCALES 

Microscopic  and  Macroscopic  Timescales 

We  treat  the  problem  by  considering  a  short  pulse  of  length  L,  where  L/c 
is  less  than  or  of  the  order  of  a  nanosecond.  Before  discussing  the  Boltzmann 
equation,  it  is  useful  to  survey  the  various  types  of  processes  involved  in 
the  problem  and  the  timescales  on  which  each  operates.  There  are  several 
timescales  relevant  to  the  dynamics: 

-  1/v,  where  v  is  an  effective  collision  frequency;  for  the  present  we 
label  all  binary  interactions  as  "collisions" .  v  may  be  less  than, 
comparable  to,  or  greater  than  o>. 


-  X/c.  This  is  the  wave  period,  the  timescale  associated  with  the 
interaction  of  a  single  electron  with  a  cycle  of  the  wave. 


L/c.  This  is  the  transit  time,  the  timescale  associated  vith  the 
interaction  of  a  single  electron  vith  the  pulse. 


®in<Vi>  V>>*  vhe«  \ 


V|vg  -  >1'  Tv2 


X/6v  ,  ,  v  and  v  ,  are 
ph  g  ph 


the  wave's  group  and  central  phase  velocities,  and  AVp^  is  the  width 
of  the  phase  velocity  spectrum  of  the  wave  packet.  These  are  the 
timescales  on  which  the  waveform  changes  due  to  dispersion.  The 
quantity  X/|v^  -  vp^l  is  the  time  required  for  a  point  of  constant 
phase  to  move  a  distance  of  order  X  in  the  rest  frame  of  the  pulse, 
i.e.  the  frame  moving  vith  the  group  velocity.  X/Av^  is  the  time 
required  for  dispersive  effects  to  produce  a  local  distortion  in  the 
shape  of  the  waveform. 

Tj  -  L^/c,  where  is  the  damping  length.  This  is  the  wave  damping  or 
attenuation  timescale. 

xa  -  h/c,  where  h  is  the  atmospheric  or  ionospheric  scale  height.  This  is 
the  timescale  on  which  the  ambient  environment  seen  by  the  wavepacket 
changes . 

Dispersion  Timescale 

For  all  cases  of  interest,  L  <<  h,  L^,  so  that  xp  <<  x^.  For  a 

narrowband  wavepacket  in  a  weakly  dispersive  medium,  we  also  have  x  «  x  . 

* 

To  obtain  a  more  quantitative  feeling  for  the  magnitude  of  the  ratio  x  /xw  we 

estimate  x  using  am  approximate  dispersion  relation.  We  note  that  for  a  wave 

^  8 

in  the  1  to  100  GHz  range,  and  assuming  an  electron  density  of  10 

3 

electrons/cm  or  less,  we  have 

(«p/«)2  «  1, 

where  wp  is  the  plasma  frequency.  In  this  case  the  dispersion  relation  is 

approximately 

2  2  2 
(kc/w)  *  1  -  («p)  /[«(»-iv)J  +  terms  of  order  (Q/<0)(wp/co)  ,  (1) 

where  2  is  the  electron  gyrofrequency,  and  v  is  the  effective  collision 

frequency.  Equation  (1)  can  be  used  to  show  that 

xwi  -  X<w2+v2)/[c(«pa>)2). 


Next  ve  estimate  x  The  width  of  the  k  spectrum  depends  on  the  shape  of  the 


"v2“ 

wave  envelope,  but  2ji/L  is  a  good  estimate  for  practical  cases.  Then  Av 


(3Vp^/3k)Ak,  and 


3v  , /3k 
ph 


ph 


can  be  calculated  from  equation  (1).  To  lowest 

2 

order  in  the  small  parameter  (w  /«)  one  finds 


v2 


L(«2+v2)/[c(«pW)2] . 


Clearly  xy^  <  x^,  since  X  <  L.  Thus  xy  =  xy^  and 

x  / x  -  (L/ X)[w  w/ (w2+v2)]2.  (2) 

p  w  p 

8 

For  a  one  nanosecond  pulse  at  1  GHz,  n  -  10  /cc,  v  *  0,  equation  (2)  yields 
-2  e 

x  /x  -  10  .  We  regard  these  as  "worst  case"  parameters:  higher  frequencies, 

p  v 

nonzero  values  of  v,  and  lover  values  of  ng  all  make  this  ratio  smaller.  In 
view  of  (2),  ve  shall  ignore  the  difference  between  the  wave  frame  (v  an<* 
the  pulse  frame  (v^)  in  the  following. 


III.  QUASISTATIC  BOLTZMANN  EQUATION 


Now  consider  the  Boltzmann  equation  in  a  frame  moving  with  the  group 
velocity  v^.  We  assume  that  the  direction  of  propagation  is  parallel  to  the  z 
axis,  and  that  spatial  gradients  in  the  x  and  y  directions  are  negligible.  In 
this  reference  frame  we  must  use  the  relativistic  Boltzmann  equation,  which  is 


3f 

3t  + 


e[E  +  (v/c)xB]-|i 


3f 

It 


coll. 


(3) 


In  the  wave  frame  the  two  dominant  causes  of  temporal  variation,  wave- 
particle  interactions  and  "collisions"  (atomic  processes)  produce  only  spatial 
variations.  In  the  absence  of  dispersion,  attenuation,  or  inhomogeneities  in 
the  ambient  environment,  ve  would  expect  steady  convective  flow  in  the  rest 
frame  of  the  pulse,  with  ambient  gas  at  all  z<0  and  a  spatially  growing 
disturbance  between  z=0  and  z=L.  We  do  not  expect  hydrodynamic  turbulence  to 
appear  within  the  pulse,  since  the  heavy  particles  cannot  respond  on  a 
nanosecond  timescale.  Moreover,  since  v  is  close  to  c,  and  x  ,  x, ,  x  are 

Z  w  u  d 

all  »  Xp,  it  follows  that 


3/3t  -  l/min(xy,  x^,  xa)  <<  v^  3/3z  -  1/x^  or  larger. 


Thus  the  quasistatic  Boltzmann  equation, 


vz 


3f 

3z 


e(  v/c)xB*-|^ 
oP 


3f 

at 


coll. 


(A) 


in  which  the  parameters  of  the  pulse,  the  ambient  medium,  and  the  distribution 
function  are  all  understood  to  be  slowly  varying  functions  of  time,  is  an 
excellent  approximation;  the  error  incurred  by  dropping  3/ 3 1  in  the  pulse 
frame  is  of  order  1%  in  the  worst  case,  as  shown  above.  Note  that  the 
electric  field  has  been  dropped  from  the  equation,  since  E  -  -  (l/c)3A/3t 
vanishes  in  the  wave  frame;  we  recover  E  when  transforming  back  to  the  lab 
frame.  In  the  wave  frame  the  total  momentum,  p,  is  a  constant  of  the  motion, 
since  we  have  only  a  static  magnetic  field.  Electrons  and  ions  enter  the 
pulse  moving  along  the  z  axis  at  nearly  the  speed  of  light;  the  difference  in 
their  velocities  is  small.  The  magnetic  field  rotates  the  large  wave  frame 
electron  velocity  vector  away  from  the  ion  velocity  vector  and  toward  the 
direction  perpendicular  to  both  z  and  the  B-field  direction  (assume  B  is 
aligned  with  the  y  axis),  and  this  appears  in  the  lab  frame  as  an  acceleration 
in  the  x-direction.  To  an  observer  in  the  lab  frame,  it  appears  as  if  the 
magnetic  force  is  weak,  and  that  the  acceleration  is  produced  by  a  strong 
electric  field  in  the  x-direction 

In  the  pulse  frame  the  "information  flow"  is  directed  from  the  head  of 
the  pulse  towards  the  tail,  for  all  quantities.  Thus,  in  addition  to  the 
reduction  of  the  dimension  of  the  system  by  one,  we  are  now  dealing  with  a  one 
point  boundary  value  problem,  and  we  can  seek  solution  algorithms  which 
involve  marching  through  the  pulse  from  head  to  tail.  In  particular,  it  is 
possible  to  develop  algorithms  involving  the  solution  of  coupled  sets  of 
ordinary  differential  equations  or  quadratures.  This  is  the  essential 
attribute  of  the  approach.  In  the  following  sections  examples  of  such 
algorithms  are  outlined.  In  principle  we  can  solve  (4),  use  f  to  calculate 
the  instantaneous  current,  electron  density,  conductivity,  etc.,  as  functions 
of  the  slowly  varying  ambient  parameters,  and  use  these  in  the  Maxwell 
equations  to  self  consistently  determine  the  slowly  varying  pulse  amplitude 
and  waveform. 
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IV.  REDUCTION  TO  SIMULTANEOUS  ORDINARY  DIFFERENTIAL  EQUATIONS 


The  lab  frame  counterpart  of  equation  (4)  can  be  written  in  terms  of  the 
coordinate  C  «  z_vpht’  v^ere  vph  =  c  *s  the  phase  velocity.  We  simply  assume 
that  the  distribution  function  in  the  lab  frame  is  a  function  of  C  and 
velocity  components  only;  justification  for  this  is  provided  by  eq.  (4).  Then 
it  follows  that,  in  the  lab  frame, 


-  V  ft  -  ><«  *  (»/c)xB(C)l-fj 


3f 

St 


coll. 


(4') 


Coulomb  or  Fokker-Plank  type  collisions  are  unimportant  on  nanosecond 
timescales,  for  n^  values  of  interest  here.  The  remaining  "collision"  effects 
can  be  expressed  as  follows: 


3f 

3t 


coll. 


J  f(C,v')  S(v, v' )  d3v' 


G(v)  f(C,v) 


The  first  term  on  the  right  describes  the  rate  at  which  particles  are 
transported  to  the  phase  space  location  (£,v)  from  all  other  phase  space 
locations  (C,v' ),  by  binary  interaction  processes.  The  second  term  gives  the 
rate  of  loss  from  (C,v)  due  to  such  processes.  The  functions  S  and  G,  which 
describe  atomic  interactions,  are  assumed  to  depend  on  the  slowly  varying 
ambient  parameters,  but  the  associated  functional  dependences  on  z,t  are 
suppressed  here,  in  keeping  with  the  quasistatic  approximation. 

Since  f  is  square  integrable,  we  can  express  it  as  a  sum  of  orthogonal 
base  functions  in  velocity  space,  with  coefficients  which  vary  with  We 
assume  that  S  and  G  can  also  be  expanded. 


f(C,v)  = 
G(v) 


S(  v,  v' )  - 


Sijk  fijk<^  *ijk(V)  =  1 


Zijk  Gijk  *ijk(v) 


^ijk^i'j'k'  Gijki'j'k'  *ijk 


Ijlj,  Sj  Jf  $j(v)  $Jr( 


J  fJ<C)  *j<v) 

J  GJ 

(v)  Vj'k'<v'> 

V') 


In  general,  the  orthogonal  functions  for  the  expansion  in  a  three 
dimensional  velocity  space  would  carry  three  indices  < *i jk.) *  s*nce  ve 

truncate  the  series  we  can  easily  rearrange  these  into  a  one  dimensional  array 
with  a  single  index  ($j)»  where  the  maximum  value  of  J  is  N,  the  total  number 


of  terras  in  the  expansion.  These  relationships  can  be  used  to  transform  (4') 
into  a  set  of  coupled  ordinary  differential  equations  for  the  coefficients 
fj(C).  In  the  simplest  case,  when  J v j  <<  c,  ve  can  neglect  the  magnetic  field 
term  in  (4')  and  set  (v^  -  v  ^)  =  -c  (these  approximations  simplify  the 
analysis,  but  neither  is  absolutely  necessary).  Substituting  the  expansion 
relations  above  into  (4')>  performing  the  velocity  integration,  re-expanding 
remaining  products  of  the  form  $j(v)#j,(v')  into  sums  of  single  *'s,  using 
recursion  relations  or  re-expansion  (depending  on  the  choice  of  orthogonal 
functions)  to  express  3*/3v  in  terms  of  *'s,  and  equating  coefficients  of  the 
orthogonal  functions  yields  the  desired  set  of  N  equations  for  the  F's: 

^  fj(t)  .  TJ,J'(Q  fJ'(0  =  0  (4") 
°r’  f(0  +  T(C)  f(0  -  0 

where  £(C)  is  the  vector  of  expansion  coefficients  of  f(C,v),  and  T  is  a 
matrix  whose  elements  are  functions  of  the  electric  field  and  the  S  and  G 
expansion  coefficients.  In  the  case  in  which  v  »  co,  the  distribution 
function  will  be  nearly  isotropic,  in  which  case  the  number  of  separate 
equations  to  be  solved  and  the  si2e  of  the  matrix  T  can  be  drastically 
reduced. 

In  any  case,  ve  solve  by  integrating  from  £  =  0  (ambient  region  ahead  of 
the  pulse)  to  C  »  L.  The  accuracy  of  the  resulting  solution  is  limited  by  the 
amount  of  computation  ve  are  willing  to  do  (manifested  by  the  number  of  terms 
retained  in  expansions,  other  approximations  employed,  etc.),  up  to  the 
accuracy  limit  of  the  quasistatic  approximation  itself  (~1X),  but  it  seems 
likely  that  good  results  can  be  achieved  with  a  reasonable  level  of  effort. 

An  alternative  approach,  in  the  same  spirit,  would  be  to  represent  the 
distribution  function  via  a  set  of  moments.  One  represents  the  distribution 
function  by  the  product  of  a  truncated  series  of  orthogonal  functions  and  a 
specified  "approximation  function"  such  as  a  Maxwellian  or  bi-Maxwellian  with 
free  parameters  (drift  velocity,  temperature,  etc.).  A  coupled  set  of 
ordinary  differential  equations  for  the  free  parameters  and  orthogonal 
function  coefficients  is  obtained  by  the  usual  method.  Convergence  of  the 
function  series  is  enhanced  when  the  "approximation  function"  is  well  chosen. 

The  coupled  ordinary  differential  equation  approach  is  likely  to  be  the 
most  practical  method  of  generating  solutions.  It  is  possible,  however,  to 


vrite  down  an  analytic  solution  of  the  quasistatic  equation  in  terms  of  a 
"scattering  operator"  F.  This  approach  is  discussed  in  the  next  section. 


V.  FORMAL  SOLUTION  OF  THE  QUASISTATIC  EQUATION 

We  note  that  in  the  absence  of  collisions  there  are  three  constants  of 

the  motion  for  a  particle  in  the  force  field  of  the  pulse,  in  the  wave  frame. 

We  have  assumed  the  wave  to  be  linearly  polarized  with  B  in  the  y-direction, 

so  that  in  addition  to  the  total  energy  the  canonical  momenta  P  and  P  are 

x  y 

solutions  of  the  homogeneous  equation  associated  with  (4).  This  means  that  ve 
can  simplify  the  operator  on  the  left  hand  side  of  (4)  by  transforming  to  a 
new  set  of  coordinates: 


where 


2’  v  V  p2  — >  z’  V  v  p 

Px  =  px  +  eAx(z)/c 


p 


y 


p 


y 


With  the  new  coordinates,  the  quasistatic  Boltzmann  equation  in  the  wave  frame 
is 


vz(z*C)  U 


const. 


3f 

3t 


coll. 


(5) 


where  C  is  a  shorthand  representation  of  the  set  of  constants  P  ,  P  ,  p.  The 

x  y 

"collision"  term  can  be  expressed  in  the  pulse  frame  as  follows: 


3f 

at 


coll. 


J  f(z,C')  S(z,C,C')  d3C' 


G(z,C)  f (z,C) 


(6) 


The  first  term  on  the  right  describes  the  rate  at  which  particles  are 
transported  to  the  phase  space  location  (z,C)  from  all  other  phase  space 
locations  (z,C'),  by  binary  interaction  processes.  The  second  term  gives  the 
rate  of  loss  from  (z,C)  due  to  such  processes.  Equations  (5)  and  (6)  can  be 
combined  and  rewritten  as  follows: 

+  g(z,C)  f (z,C)  =  I(z,C)  (7) 
=  Jf(z,C')  s(z,C,C' )  d3C'  (8) 


where 


df(z,C) 


I(z,C) 
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g(z,C) 


(9) 


,-l 


(vz(z,C)i  G(z , C) 


,-l 


s(z,C,C' )  =  [vz(z,C)]  S(z,C,C' ) 


Note  that  3/3z  has  been  replaced  by  d/dz.  Since  the  C's  are  to  be  held 
constant  under  the  dif ferentation  by  z,  ve  can  treat  (7)  as  an  ordinary  first 
order  differential  equation  in  z,  if  ve  regard  I(z,C)  as  a  known  function.  In 
terms  of  I(z,C),  the  solution  to  (7)  is 


f (z,C)  =  e-Jg(z' .C)dz'  [  f(Q)C)  +  jz  I(z,fC)  eJg(z",C)dz"  dz, 


(10) 


Nov  define  the  function  space  operator  T  such  that 


rf  =  Jq  ;  f(z',C')  s(z' ,C,C' )  d3C'  dz' 


(ID 


Eq.  (11)  becomes 

e/g(z' >C)dz'f(2>C)  =  f(0|C)  +  re/g(z' ,C)dz'f 


(12) 


Ve  can  iterate  (12)  to  obtain 


!Jg(zSC)dz'f(2jC)  a  [1  +  r+  r2  +  r3  +  ...j  f(Q,C) 


(13) 


We  can  formally  "sum"  (13)  and  express  the  solution  to  (12)  in  the  form 


f(z,C)  *=  e-Jg(z',C)dz'(1  _  j.j-1 


f(0,C) 


(14) 


which  can  also  be  obtained  directly  from  (12). 

The  operator  T  can  be  interpreted  as  a  "scattering  operator".  In  the 
absence  of  "collisions",  the  exact  solution  of  the  quasistatic  equation  is 
given  by  the  first  terra  of  (13);  the  vavefield  effects  are  carried  in  the 
constants  of  motion  C(z,p).  In  a  veak  scattering  situation,  where  electrons 
do  not  undergo  more  than  one  atomic  interaction  during  their  transit,  T 
operating  on  the  primary  flux  f(0,C)  describes  the  contribution  to  f  from 
electrons  scattered  into  the  interval  (C,C+dC)  after  a  single  interaction,  and 
ve  expect  that  (l+T)f(0,C)  vill  be  an  adequate  solution.  Similarly,  ve  can 
interpret  I^ftOjC)  as  the  contribution  from  electrons  that  enter  (C,C+dC) 
after  being  scattered  tvice.  rnf(0,C)  is  interpreted  in  analogous  fashion. 


The  effects  of  losses  from  the  (C,C+dC)  interval  are  described  in  closed  form 
by  the  e-^2  term.  In  general,  if  these  interpretations  are  correct  and  if 
L<.  is  the  mean  free  path  for  atomic  processes  which  add  particle-  to  the  phase 
space  interval,  then  ve  expect  that  the  series  in  (13)  has  converged  after  n 
terms  if  n  is  much  larger  than  the  expectation  value  fr r  the  number  of 
scatterings  experienced  by  an  electron  in  moving  from  0  to  z,  i.e.  n  >>  z/L,,. 
To  illustrate  and  test  these  ideas  consider  a  simple  model  problem  which  can 
be  solved  exactly.  Let  f(z,C)  be  a  function  of  z  only,  f  -  F(z),  inside  a 
region  R  of  C-space  with  volume  V^,  and  zero  outside  this  volume.  Also  take 
g(z,C)  -  1/LG  and  s(z,C,C')  -  l/(Lg  V^).  Equation  (7)  then  has  the  solution 


f(z,C)  -  F(z)  ^  F(0)  exp[z/Lg  -z/LG]  *=  F(0)  exp[z/L']  for  C  in  R 

f(z,C)  =  0  for  C  outside  R 


where  L'  «  LgLG/(LG  -  Lg) 

If  Lg  <  LG,  then  L'  is  positive.  Ve  have  an  avalanche  if  0  <  L'  <<  L. 

Now  let  us  consider  (13).  For  the  model  problem,  we  find 

f1  F(0)  -  (z/Ls)n/n!  F(0)  (15) 


With  (15),  (13)  is  easily  summed  to  yield  the  correct  solution  to  the  model 
problem.  It  follows  from  (15)  that  the  ratio  of  the  adjacent  terms  rn/rn-'*'  in 
(13)  is  (z/Ls)/n,  so  that  the  convergence  condition  is  the  same  as  that 
inferred  above.  Note  that  the  series  always  converges,  even  in  the  case  of  an 
avalanche. 

Ve  can  express  the  scattering  operator  as  a  matrix  by  expanding  in 
orthogonal  functions,  as  before.  Ve  assume  the  orthogonality  condition  has 
the  form  (a  slight  modification  of  the  analysis  is  needed  if  there  is  a  weight 
function) 


I  d3C 


Eq.  (12)  can  be  written  in  the  form 


Q(z,C)  =  f(0,C)  +  FQ(z,C) 


(16) 


where 


Q(z,C)  =  e 


Jg(z' ,C)dz' 


f(z,C) 


Differentiation  of  (16)  with  respect  to  z  yields 


-  J  s(z,C,C' )  Q  d3C' 


Now  expand  s  and  Q  in  orthogonal  functions,  and  use  the  orthogonality  relation 
to  eliminate  the  integral.  We  get 

dQ/dz  -  s  Q  =  0  (17) 


where 


Qj(2)  -  S  Q(z,C)  *j(C)  dJC 

Sj  j,(*)  -  //  s(z,C,C' )  *j(C)  *j,(C')  d3C  d3C' 


The  solution  to  (17)  is 


CO 

Q(z)  »  exp[ Js(z' )dz' ]  Q(0)  *  f  £  [ Js(z' )dz' ]n/n! ]  )  Q(0) 

t  o  -  1 

or,  since  $*Q  *  exp[ Jg(z' ,C)dz' ]  f(z,C)  and  Q(0)  *  f(0), 


f(z,C)  ■  4-exp[  Js(z'  )dz'-IJ‘g(z'  ,C)dz'  ]  •  f (0) 


Eq.  (18)  is  similar  in  form  to  the  solution  of  the  model  problem  considered 
above.  For  computation,  one  would  replace  the  matrix  part  of  the  exponential 
in  (18)  by  its  series  expansion  ((18)  is  then  the  matrix  version  of  (13)]  and 
transform  back  to  the  lab  frame.  In  its  present  form,  (18)  is  more  cumbersome 
for  computational  purposes  than  (4"),  especially  in  the  collisional  regime, 
but  it  is  possible  that  further  analysis  will  lead  to  approximate  algorithms 
useful  in  the  weakly  collisional  regime  where  the  series  expansion  converges 
rapidly. 

In  the  next  section  we  develop  in  some  detail  the  equations  needed  to 
implement  a  version  of  the  coupled  ordinary  differential  equation  approach 
discussed  in  section  IV. 
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VI.  COUPLED  ORDINARY  DIFFERENTIAL  EQUATION  CALCULATIONS 


The  General  Formalism 

In  the  laboratory  frame,  the  Boltzmann  equation  can  be  written  as 


<vz  -  vH  '  » [E(t)  *  v/c  ’‘B<0' 


3f 

3v 


3f 

at 


coll. 


(19) 


where  we  have  used  m  »  electron  mass,  £  * 


z-v  ,  t  and  3f/3t 
ph 


In 

this  section,  we  solve  the  above  equation  by  reducing  it  to  a  system  of 
ordinary  differential  equations.  For  the  purpose  of  this  report,  which  is  to 
demonstrate  the  formalism,  we  will  use  a  simplified  scattering  model. 
However,  the  basic  technique  is  not  limited  to  any  particular  scattering 
model. 

we  neglect  the  magnetic  field  term  and  set 


(v 


Since  v  ^  »  c  and  |v|  <<  c, 


vphJ 


-c  in  the  above  equation. 


Ve  also  assume  that  the  system  is 
nonrelativistic.  The  coordinate  system  is  such  that  z  is  the  propagation 
direction  of  the  microwave  pulse  and  the  electric  field  is  in  the  x  direction. 
Then,  the  approximate  equation  can  be  written  as 


3f 

+ 


inc  Ex«> 
e 


1  3f 
c  3t 


=  0 


(20) 


The  last  term  describes  the  time  rate  of  change  of  the  distribution  function 
due  to  collisions.  In  the  present  problem,  we  have  partially-ionized  plasmas 
so  that  the  short-ranged  collisions  of  electrons  with  neutral  particles 
dominate.  Coulomb  collisions  are  unimportant.  Thus,  the  collision  term  takes 
on  the  general  form 


-  JdVf(C,  v'  )S(v,  v')  -  G(v)f(C,v) ,  (21) 

where  the  first  term  gives  the  time  rate  of  increase  in  the  distribution 
function  f(C,v)  due  to  collisions  at  all  v'  and  the  second  term  represents  the 
time  rate  of  loss  of  electrons  at  (C,v)  due  to  collisions  taking  place  at 
(C,v). 


3f 

3t 


12 


The  plasma  distribution  function  f  is  taken  to  depend  on  £  and  the 
spherical  velocity  coordinates  v  and  6,  where  9  is  the  polar  angle.  Ve  define 
v  ■  v^/v  -  cos9  and  use  a  separable  form  of  f  and  the  following  set  of  basis 
functions: 


and 


♦j(»)  •  (^-)1/2  v*1sin  j  - 


Va)  *  prrS11'2  V">'  »  -  2'3> 


The  normalization  relations  are 


J  dv  v2^<v)(j^(v)  -  5jk 


and 


j!, 


dvP  (v)P  (v)  .  Sm, 
m  n  ZJl  mn 


Here,  vq  is  the  maximum  electron  velocity,  and  Pm  (w)  is  the  usual  Legendre 
polynomials.  Ve  then  exand  f,  S  and  G  as  follows: 


Then,  (20)  can  be  written  as 
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i ,  1  j ,  m 


where  -  d+./dv  and  P^'(w)  »  dP^/dw.  If  we  multiply  all  terms  by  the  basis 
functions  ♦j(v)P]n(w)  and  integrate  over  velocity  space,  then  it  is  clear  that 
a  system  of  ordinary  differential  equations  results. 

A  Simple  Collision  Model 

In  order  to  solve  (22),  we  need  to  specify  the  collision  terms.  While 
this  equation  can  be  used  for  any  scattering  cross-sections,  we  will 
illustrate  the  method  and  its  capabilities  by  using  a  simple  collision  model 
and  solving  the  equation  directly.  In  this  paper,  we  assume  that,  after  a 
collision  with  a  neutral  particle,  an  electron  has  an  equal  probability  of 
being  scattered  into  any  solid  angle.  Then,  we  have 


G(v)  -  n.S 

a  s  s 


S(v,v')  *  if  v_1E  o  (V)  &  v'  -  V  +e  , 


m 


TO 


WKW5 


cSw*! 
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vhere  s  is  the  index  for  various  collision  processes,  n  is  the  ambient 

neutral  density  and  is  the  scattering  cross-section  (integrated  over  the 

solid  angle)  for  the  collision  process  s.  The  quantity  eg  is  the  normalized 

threshold  energy  for  the  interaction  s;  an  electron  vith  incoming  energy  of 
2 

(m@v  /2)  loses  (mgcg/2)  as  a  result  of  the  collision.  For  each  interaction  s, 

if  the  incident  energy  is  less  than  s  ,  then  no  scattering  takes  place  so  that 

2  s 
a  (v)  -  0  for  v  <  e  .  We  define 
s  s 

<*T(v)  *  (v)  (25) 

s 

Using  (23)  and  (24),  it  is  straightforward  to  calculate  the  expansion 
coefficients  and  S^jm.  By  multiplying  all  the  terms  in  (22)  by 

♦j(v)Pm(v)  and  integrating  over  the  velocity  space,  (22)  can  be  reduced  to  the 
following  system  of  differential  equations, 


*  *  rh*uK -i<«] 

*  v°  ■  °- 

where  i,j,k  *  1,2,3,..  and  l,m,n  =  0, 1,2,3 .  The  collision  term 

Xjm  can  be  written  as 

V«  ■  (r)(rko£fio<^i  - 

0  i  o  IX. 


In  the  above  expressions,  the  following  definitions  have  been  used: 


A..  *  [  dx  x~^sin(ix)sin(jx), 
ij  J0 


L.  .  ■  i  f  dx  cos(ix)  sin(jx)  -  ij  [l  -  (-l)^+i](j2  -  i2)-1 
o  J0  L  J 
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dx  x 


sin(ix)sin(kx)sin( jx) . 


then  we  find 
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8n  it  .1/2 
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n  v  1  /  2 
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using  the  definitions  of  3..  and  S.,  it  can  be  shown  rigorously  that 


In  (27),  Xjm  (O  describes  the  changes  in  the  distribution  function  due  to 
the  creation  of  ionization  secondary  electrons.  This  term  will  be  calculated 
later  according  to  a  simplified  model.  Note  that,  because  of  the  isotropic 
scattering  assumption,  ve  have 


S . ,  .  -  (2-  1  n  8,  8  3 .  . 

llim  (v  J  a  lo  mo  ii 


Gil-<4")  rj  "a5lo2i 
o 

vhere  8  denotes  the  8-function.  For  more  realistic  scattering  cross-sections 

of  the  form  cr  (v,0),  the  summations  in  (26)  would  be  more  complicated.  Note 

that  conservation  of  particles  requires  that  JdJv(3f/3t)|  ,,  »  0  in  the 

( 1  \  c  • 

absence  of  ionization  (Xjffl  *  0)  for  any  distribution  function  f.  If  we 

define 


K*  f  If 


2v  ® 
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z  ^  \ 

j.l  J 


for  any  choice  of  <*s(v),  excluding  ionization.  Thus,  T  ■  0  exactly  for  any  f 
and  any  scattering  processes. 

Solution 

We  have  solved  the  set  of  differential  equations  (26)  by  numerical 
integration.  For  the  present  work,  ve  used  21  scattering  processes  to 
calculate  3^  and  5^:  they  are  ionization,  vibrational  excitations  of  N ^ 
molecules,  electronic  excitations  of  Nj  molecules,  and  momentum  transfer.  For 
the  purpose  of  numerical  computation,  it  can  be  shown  that 


vhere 


\j  *  J  dv  vfft(v)sin(  J^vjsingiv). 


Using  (29),  ve  can  rewrite  (27)  as 


V«  ‘  £i„<«Aji] 


(30) 


This  expression  eliminates  the  double  summation,  providing  better 
accuracy.  For  the  purpose  of  illustrating  the  procedures  involved  in 
calculating  the  contribution  of  the  secondary  electrons  created  by  ionization, 
ve  assume  a  simple  model  in  which  the  secondary  electrons  are  produced  with 
zero  energy.  The  result  is  expected  to  depend  on  the  distribution  of  the 
secondary  electrons,  and  for  serious  calculations  one  would  use  a  more 
realistic  distribution.  Pitchford^,  for  example,  has  assumed  that  the 
incident  and  secondary  electrons  share  the  energy  equally  following.  In  this 
paper,  we  will  not  discuss  the  detailed  consequences  of  the  model  dependence; 
a  task  for  future  studies.  If  ve  assume  that  all  secondary  electrons  are 
created  with  zero  initial  energy,  then 


3f 

at 


(i) 

coll. 


S3(v)  J  G(i)(v')  f(C,  v')  dV, 


where  the  superscript  (i)  means  that  only  the  ionization  cross-section  is 
used.  This  expression  gives  the  time  rate  of  change  of  the  distribution 
function  due  to  the  ionization  secondary  electrons.  Using  this  expression,  we 
find 


(2  Jin 


(cv 


(31) 


where  G .  is  obtained  from  the  definition  of  G.  following  (27)  by  replacing 

1  ( i)  1 

(Tj(v)  with  a  (v),  the  ionization  cross-section. 

Our  code  solves  the  set  of  equations  represented  by  (26),  with  the 
collision  terms  Xj  m ( C)  given  by  (30)  and  (31)  or  other  expressions 
representing  different  assumptions  about  the  nature  of  the  "collisions".  In 
particular,  the  assumption  that  all  the  secondary  electrons  are  created  with 
zero  energy,  discussed  above  as  a  simple  illustration  of  the  method,  is 
numerically  troublesome  and  somewhat  unphysical:  in  practice  the  newly  created 
electron  population  must  have  a  finite  temperature.  The  electron  distribution 


function  can  be  reconstructed  and  displayed  in  phase  space,  or  moments, 
conductivities,  and  other  quantities  of  interest  can  be  calculated  directly 
from  the  function  space  representation. 

For  the  purpose  of  evaluating  the  quantities  in  (26),  ve  define  equally 
spaced  grids  for  the  velocity  and  angular  coordinates.  Then,  all  quantities 
including  the  cross-sections,  and  are  defined  and  computed  on 
these  grids.  For  diagnostic  purposes,  ve  check  particle  and  energy 
conservation. 

Figures  1,  2,  and  3  illustrate  the  type  of  results  one  obtains  with  this 

method,  as  applied  to  a  sample  problem:  E  =  100  statvolts/cra. ,  w/2rt  =  30 

wave  2 

GHz,  and  ve  have  an  atmosphere  of  Nj  vith  a  density  of  10  /cm  .  Secondary 
electrons  are  created  vith  a  thermal  spread  of  2  ev,  in  this  test  calculation. 
Figure  1  represents  the  drift  of  an  initially  Maxvellian  distribution  function 
vith  20  ev  temperature  under  the  action  of  a  vave  electric  field  only  (the 
vave  amplitude  has  been  reduced  to  25  statvolts/cm.  to  avoid  runaway):  atomic 
processes  are  turned  off.  In  Figure  2  ve  see  the  evolution  of  the  same 
initial  distribution  function  in  the  absence  of  a  vave  electric  field,  but 
vith  atomic  processes  operative.  A  tabulation  of  the  cross  section  data  used 
is  provided  in  the  Appendix.  The  "ledge"  seen  to  form  on  the  shoulders  of  the 
Maxvellian  is  due  to  vibrational  excitations  vith  cross  sections  sharply 
peaked  near  2.5  ev.  The  groving  central  peak  is  primarily  due  to  the  cooling 
effect  of  the  various  excitation  processes;  i.e.,  these  are  mostly  electrons 
scattered  dovn  from  higher  energies.  Ionization  of  neutrals  by  hot  tail 
particles  also  contributes  to  the  growth  of  the  central  peak.  In  Figure  3, 
the  vave  (100  statvolts/cm.  amplitude)  and  the  atomic  processes  act 
simultaneously,  and  one  can  observe  drift  motions,  the  buildup  of  electron 
density  due  to  ionization,  etc.  Note  that  the  distribution  function  quickly 
becomes  non-Maxvellian. 

VII.  SUMMARY  AND  COMMENTS 

Ve  have  discussed  an  approach  to  the  solution  of  the  Boltzmann  equation 
suitable  for  the  study  of  short  intense  microwave  pulses  interacting  vith  a 
background  atmosphere  in  various  parameter  regimes,  and  a  specific  example 
appropriate  in  the  collisional  regime.  The  coupled  ordinary  differential 
equation  approach  appears  to  be  viable.  Continuing  studies  and  comparisons  of 
this  approach  vith  others  are  expected  to  expose  the  advantages  and 
disadvantages  of  the  method,  and  lead  to  refinements. 
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Figure  1.  This  Figure  shows  the  drift  of  an  initially  Maxwellian  distribution 

function  with  20  ev  temperature  under  the  action  of  a  wave  electric 

field  (E  «  25  statvolts/cm. ,  w/2n  »  30  GHz)  only:  atomic 

wave 

processes  are  turned  off.  The  four  subfigures  represent 

"snapshots"  of  the  distribution  at  different  wave  phases,  as  the 
wave  pulse  moves  past  an  element  of  the  electron  gas.  The  phases 
displayed  are:  (top  left)  =  X/4,  (top  right)  ^  *  X/2,  (bottom 

left)  C2  -  3X/4,  (bottom  right)  =  X,  where  X  is  the  wavelength. 


vith  cross  sections  sharply  peaked  near  2.5  ev.  Distributions  are 
displayed  at  "times"  C^/c,  i  =  1,2, 3, 4. 
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APPENDIX:  CROSS  SECTION  TABLES 

Cross  section  data  used  in  the  sample  calculation  described  in  the  text 
is  tabulated  belov.  The  threshold  energy  for  each  process  is  given  in 
parenthesis. 


N2  Rotational  Cross  Section 


(Eq  -  .02  ev) 


Electron  Energy  (eV) 


7.2 ( — 18 ) 
2.1(-17) 

2  -  7 ( —17 ) 
4.8(-17) 

1 ♦ 7 ( —16) 
6.2(— 16) 
6.5(-16) 
6. 1 ( — 16 ) 
8. 28(-16) 
7.2(-16) 
6.1 (—16 ) 
2.7 ( — 16 ) 
4.5(-17) 

0 
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N2  Vibrational  Cross  Section  (v  -  1) 

ff1  (E1  *  0.29  ev) 

Electron  Energy  (eV) 

o(cm^) 

0.3 

1 . 32(— 19) 

0.75 

6.4(-19) 

1.44 

1 . 22(— 17 ) 

1.65 

1.84(-17) 

1.74 

2.30(-l7) 

1.85 

4.3(-17) 

1.92 

1 . 01 (—16 ) 

1.98 

2 . 60(-16) 

2.02 

3.0(-16) 

2.09 

2.1(-16) 

2.12 

1.46(-16) 

2.18 

7.8(-17) 

2.24 

1 . 38 ( — 16 ) 

2.27 

2  -  36 (  — 16 ) 

2.31 

2.8(-16) 

2.38 

2 . 18( —16 ) 

2.41 

1. 39(-16) 

2.47 

1.07(-16) 

2.54 

1 . 2 ( —16 ) 

2.57 

1.78(-16) 

2.61 

2.0(-16) 

2.66 

1 . 66 ( — 16 ) 

2.71 

3.6(-17) 

2.73 

8. 8(-17) 

2.76 

5. 6(-17) 

2.79 

7 .0(-17) 

2.83 

3  -  68 (  —  17) 

3.07 

7 . 68(-l7) 

3.27 

3.98(-17) 

3.60 

3. 38(-17) 

4.0 

0 
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Nj  Vibrational  Cross  Section  (v  ■  3) 


<j3  (E3  -  0.87eV) 
Electron  Energy  (eV) 


g(cm^) 


1.80 

0.0 

1.87 

3.0(-18) 

1.97 

4.0(-17) 

2.02 

7  -  8 ( — 17 ) 

2.05 

1 . 17 ( -16 ) 

2.12 

1 . 39( — 16 ) 

2.17 

1 . 72 ( — 16 ) 

2.22 

1.76<— 16) 

2.27 

1.38(-16) 

2.32 

7  -  2 ( —17 ) 

2.36 

2.24(— 17) 

2.40 

3.0(-18) 

2.47 

4 . 2 ( — 17 

2.52 

6.9(-17) 

2.57 

8.7(-17) 

2.6 

9.6(-17) 

2.67 

6.9(-17) 

2.72 

2. 7 ( —18 ) 

2.77 

3 . 0( — 18 ) 

2.82 

3.0(-17) 

2.92 

5.1 (—17 ) 

3.02 

2.1(-17) 

3.13 

1.5(— 17) 

3.21 

1.94(-17) 

3.32 

3. 0(-18) 

3.35 

0.0 

27 


■ 


■V\ 


mid 


W-'J 


N2  Vibrational  Cross  Section  (v  -  4) 


<t4  (E4  -1.16  eV) 


Electron  Energy  (eV) 

u(cmZ) 

2.0 

0.0 

2.07 

3 . 0  ( — 18 ) 

2.13 

6.0(-17) 

2.19 

7 . 94 (—17) 

2.20 

1.19(-16) 

2.28 

1.49(-16) 

2.32 

1 . 27 ( — 16 ) 

2.37 

9.74(-17) 

2.4 

8. 1(— 17) 

2.5 

3.0(-18) 

2.56 

3.0(-18) 

2.62 

3.0(-17) 

2.67 

4 . 5 ( — 17 ) 

2.71 

7.8(-17) 

2.77 

5 . 7 ( —1 7 ) 

2.82 

2.7<— 17) 

2.86 

3.0(-18) 

2.94 

1.2(-17) 

3.0 

2.4(— 17) 

3.15 

1.8(-17) 

3.21 

3.0(-18) 

3.25 

0.0 

Nj  Vibrational  Cross 

Section  (v  -  5) 

<?5  (e5  -  1. 

45  eV) 

Electron  Energy  (eV) 

2.0 

0.0 

2.13 

3.0(-18) 

2.17 

3.9(-17) 

2.22 

5.7(— 17) 

2.32 

8.1(-17) 

2.38 

9. 3(-17) 

2.42 

9.9(— 17) 

2.48 

9.9(— 17) 

2.52 

8.1(1-17) 

2.57 

4.2(— 17) 

2.62 

3.0(-18) 

2.68 

3.0(-18) 

2.72 

1.5(-17) 

2.82 

4.2(-17) 

2.92 

4 . 2  ( — 1 7 ) 

3.0 

3.0(-18) 

3.12 

1.5(-17) 

3.22 

2.24(-17) 

3.32 

3.0(— 18) 

Nj  Vibrational  Cross  Section  (v  *  6) 


a6  (Eg  *  1.74  eV) 


Electron  Eners 


o-(cm^) 


3.0(— 18) 

4.2(-17) 

7.0(-17) 

9.0(-17) 

9.06<-17) 

9 . 06 (—17) 

8.44(-17) 

4 . 22 (  —  17) 

1.81(-17) 

1-5 ( — 17 ) 

2.7(-17) 

1  -  2 ( — 17 ) 

0.0 


Vibrational  Cross  Section  (v  =  7) 
ffy  (Ey  -  2.03  ev) 


Electron  Energy  (eV) 


g(cm^) 

0.0 

3.0(-18) 

1.66(-17) 

2.4(-17) 

3 . 02(— 17 ) 
3.92(-17) 
5. 12(-17) 
6 . 04 (—17) 
5. 12 ( — 17 ) 
3.92(-17) 
3.0(-17) 
1.2(— 17) 
3.0(-18) 
1.2(-17) 
l-5(— 17) 
9.0(-18) 
0.0 


N2  Vibrational  Cross  Section  (v  -  8) 


ff8  (Eg  =2.32  eV) 

2 

Electron  Energy  (eV)  g(csn  ) 


2.50 

0.0 

2.57 

3.0(-18) 

2.62 

1 . 8 ( — 17 ) 

2.66 

2.42(-17) 

2.71 

2.42(-17) 

2.78 

1.96(-17) 

2.84 

2.4(-17) 

3.0 

9.0(-18) 

3.04 

3.0(-18) 

3.15 

3.0(-18) 

3.3 

9.0(-18) 

3.4 

3.0C-18) 

3.5 

0.0 

Excitation  Cross  section 

For  N2(A3Z) 

<Jg  (Eg  =  6.17 

eV) 

Electron  Energy  (eV) 

a  (cm^) 

6.17 

0.0 

6.50 

1 . 28 ( — 18 ) 

7.0 

1.88(-18) 

8.0 

6.3(-18) 

9.0 

1 . 21(-17) 

10 

1 . 68( — 17 ) 

11 

2.22(-17) 

12 

1.92(-17) 

13 

1 • 5 ( — 17 ) 

14 

1-2 ( — 17 ) 

15 

1 . 0( — 17) 

16 

8.8(— 18) 

18 

6 . 7 ( — 18 ) 

20 

5 . 4 ( — 18 ) 

25 

3.8(— 18) 

30 

2 . 5 ( — 18 ) 

35 

1.7(— 18) 

40 

1.3 ( — 18 ) 

Above  40  eV  use 

8.06  x  10-14 

a  «  - b - 
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Excitation  Cross  section  For  ^(B  H) 

*10  <E10  "  7‘35  eV) 

2 

Electron  Energy  (eV)  g(cm  ) 


7.35 

0.0 

7.5 

4.0(-19) 

8.0 

5.4(-18) 

9.0 

1 . 4  ( —  1 7 ) 

10.0 

2.25(— 17) 

11.0 

2 . 78( —17 ) 

12.0 

2.99(-17) 

13.0 

2 . 97 ( — 17 ) 

14 

2.71 (—17 ) 

15 

2 . 41 ( — 17 ) 

16 

2. 16(-17) 

17 

1 . 95(— 17 ) 

18 

1 . 79(— 17 ) 

19 

1 . 66(-17) 

20 

1 . 56(-17) 

22 

1 . 42 (—17) 

24 

1 . 30 (—17) 

26 

1.20(-17) 

28 

1 . 10( — 17 ) 

30 

1 . 01 ( — 17 ) 

32 

9 . 2(-18) 

34 

8.4(-18) 

40 

5 . 1(— 13) 

Above  40  eV  use 

<t  =»  3.3  x  10~13  ~ 

E 
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Excitation  Cross  section  For  ^(v  A) 

*ii  (En  “  7,36  eV) 

Electron  Energy  (eV) 


g(cm3) 


7.36 

0.0 

8 

2 . 7(-18) 

9 

7 . 4(-18) 

10 

1 . 2  ( —  1 7 

11 

1.66(-17) 

12 

2 . 13 ( — 1 7 ) 

13 

2 . 6(— 17 ) 

14 

3.06(-17) 

15 

3 . 51 ( —17 ) 

16 

3.8(— 17) 

17 

3. 76(-17) 

18 

3 . 5(— 17 ) 

19 

3 . 09(— 17) 

20 

2 . 65 ( — 17 ) 

22 

1.97(— 17) 

24 

1 -53<— 17) 

26 

1.26(-17) 

28 

1.08(-17) 

32 

8 . 35 ( — 18 ) 

36 

6.6(-18) 

38 

5.9(-18) 

40 

5 . 2( — 18 ) 

50 

3 . 0( —18 ) 

Above  50  eV  use 
i 


a  -  3.75  x  10"13  K 


Excitation  Cross  section  For  NjCB'  E) 
<*12  (Eu  =  8.16  eV) 


Electron  Energy  (eV) 


2 

q(cm  ) 


Above  50  eV  use 

q  -  2.37  x  10'13  “t 
EJ 


1 . 6(-18) 
3.5(-18) 

5 . 5  ( - 18 ) 

7 . 4  ( — 18 ) 
9.4(-18) 
1. 13 ( - 1 7 ) 
1 . 2  5  (  —  1 7 ) 
1 . 14(-17) 
9 . 2 ( — 18 ) 

7 . 3(-18) 
6. 1 ( — 18 ) 

5 . 4(-18) 
4. 7(-18) 

4 . 3(-18) 
3. 9(-18) 

3 . 4(-18) 
2.9(-18) 
2.4(-18) 

2 . 2(-18) 

1 . 9(-18) 


Excitation  Cross  section  For  N2(a'  i) 


a13  ^E13  “  8-4  eV^ 
Electron  Energy  (eV) 


g(cm^) 


8.39 

0.0 

9.0 

1 . 0( — 18 ) 

10 

2.7(— 18) 

11 

4. 5(-18) 

12 

6.2 ( — 18 ) 

13 

8. 0(-18) 

14 

9.6 ( — 18 ) 

15 

1.04(-17) 

16 

8.5(-18) 

17 

6.4 ( — 18 ) 

18 

5.2(-18) 

19 

4.5(-18) 

20 

4. 1(-18) 

22 

3 . 4(-18) 

26 

2.7(— 18) 

30 

2.3 ( — 18 ) 

34 

2.0(-18) 

38 

1 . 87 ( -18) 

42 

1.84(-18) 

46 

1.82(-18) 

50 

1.80(-18) 

Beyond  50  eV  use 
.  -  9  «  1017 (P 


Excitation  Cross  section  For  ^(a  B) 
»14<e14  =  8.55  eV> 

2 

Electron  Energy  (eV)  g(cm  ) 


8.54 

0.0 

9 

9.5(-19) 

10 

3 . 0(— 18 ) 

11 

4. 95(-18) 

12 

7 ,0(-18) 

13 

9 . 0( —18) 

14 

l.K-17) 

15 

1.3(-17) 

16 

1 . 43 ( — 17 ) 

17 

1.50(-17) 

18 

1.49(-17) 

19 

1 . 43(— 17) 

20 

1 . 38 ( —17) 

22 

1 . 26 (  —  17) 

26 

1 . 14( —17 ) 

30 

1 . 02 ( — 17 ) 

34 

9. 2(-18) 

38 

8 . 8(— 18) 

42 

7 . 6(-18) 

46 

6 . 9(— 18 ) 

50 

6. 3(-18) 

Beyond  50  eV  use 
3  x  1016 


Excitation  Cross  section  For  Njtv  A) 
*15  (E15  =  8,87  eV) 

2 

Electron  Energy  (eV)  g(cm  ) 


8.89 

0.0 

9 

2 . 0(-19) 

10 

3 . 7 ( —18 ) 

11 

7 • 1 ( — 18 ) 

12 

9. 9(-18) 

13 

1 . 77 (-17 ) 

14 

1 . 1 5  ( —  1 7 ) 

15 

1 - 0(— 17 ) 

16 

8. 1 ( - 18 ) 

17 

6 . 6 ( —18 ) 

18 

5. 6(-18) 

20 

4. 3(-18) 

24 

3 . 2 ( — 18 ) 

28 

2 . 6(-18) 

32 

2. 1 ( —18 ) 

36 

1 . 6( — 18 ) 

40 

1 . 3 ( — 18 ) 

44 

1 . 0( —18 ) 

48 

8 . 0(-19) 

50 

7.0(— 19) 

Beyond  50  eV  use 
1.87  x  IQ-15 


'■y< 


i  ikj  iSaiVu 
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Nj  Dissociation  Cross  Section 


ff16  (E16  = 

Electron  Energy  (eV) 


10.0  eV) 

g( cm2 ) 


9.7 

0.0 

10.0 

1 . 2  7 ( — 18 ) 

10.5 

1.87(-18) 

11 

2 . 49 ( — 18 ) 

12 

4 . 92 ( — 18 ) 

13 

1 . 1 2  ( —  1 7 ) 

14 

2. 43 (—17) 

15 

3.58(— 17) 

16 

4.75(-17) 

18 

6.82(— 17) 

20 

8 . 56 ( — 17 ) 

22 

1.05(-16) 

26 

1 . 36 ( — 16 ) 

30 

1 . 58(-16) 

35 

1.67(— 16) 

40 

1.73(— 16) 

50 

1 . 75(— 16) 

60 

1. 79(— 16) 

70 

1 . 80(-16) 

90 

1 . 71 ( —1 6 ) 

100 

1.64(-16) 

After  100  eV  use 


a 


1.64  x  1Q~14 
E2 


V  J 


/V  V  v  ■  ■ 


11.03 

0.0 

11.5 

8 . 0( — 18 ) 

12 

1 ,46(-17) 

13 

2  -  98 (  —  17) 

14 

4.43(-17) 

15 

3 . 89 ( — 17 ) 

16 

2.84(-17) 

17 

2.34(-17) 

18 

2.02(-17) 

19 

1.81(-17) 

20 

1 . 65(-17) 

22 

1 . 39(— 17 ) 

24 

1 . 18  < — 17 ) 

28 

8 . 6(-18) 

32 

6 . 6(-18) 

36 

5 . 2 (  — 18 ) 

40 

4 . 2 ( —18 ) 

44 

3.4(-18) 

48 

2.8(-18) 

50 

2 . 6( — 18 ) 

Beyond  50  eV  use 

1.42  x  10~1A 
*"  E2’2 


.  f  *  1 *  .  1  •  -  ,  '  .  *  ,  ’  »*  *  <  «  *  *  ' 

■'  J>r 


r  _  "  J 


■» *.  '« ' »_ • . b ** 


•  -  \*  -j»  \ 


«W 


IonizatioD  Cross  Section 


a19  ^E19  “  15. 6 
Electron  Energy  (eV) 

16 

16.5 

17 

17.5 

18 

18.5 

19 

19.5 

20 

20.5 
21 
22 

23 

24 

25 

26 
28 
30 
32 
36 
40 
45 
50 
60 
70 
85 
100 
120 
160 
180 
200 
300 
400 
500 
600 
800 
1000 


Beyond  1000  eV  use 

<r  «  2~-°'9— f  10 -  [log  (0.08E)  -  log  (1  - 


ev) 

2 

g(cm  ) 

2.1 (—18 ) 
4.65(-18) 
7 . 12(-18) 
9.84(-18) 
1.29(-17) 
1 . 63(— 17 ) 
1 . 98 ( — 17 ) 
2. 3(-17) 
3.07(-17) 
3.07(-17) 
3 . 45 (—17) 
4. 17(-17) 
4.91(-17) 
5.65(-17) 
6. 38(-17) 
7. 13(-17) 
8. 74(-17) 
1.03(-16) 
1.15(-16) 
1 . 38(-16) 
1.57(-16) 
1. 77(-16) 
1 . 93(— 16) 
2 . 18(-16) 
2 . 33(-16) 
2.46(-16) 
2.52(-16) 
2-52(— 16) 
2.42(-16) 
2.35(-16) 
2.27(-16) 
1.92(— 16) 
1 . 66 ( — 16 ) 
1.45(-16) 
1. 29(-16) 
1.06(-16) 
9 . 22 ( - 17 ) 
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Moaentua  Transfer  Cross  Section 


m 


Electron  Energ 

0.01 

0.014 

0.022 

0.032 

0.065 

0.1 

0.23 

0.33 

1.0 

1.5 

1.6 
1.8 
2.0 
2.2 
2.6 
2.8 
3.0 
3.3 
3.6 
4.0 
4.5 
5.0 
6.0 
7.0 
9.0 
10 
15 
20 
30 
50 
100 
200 
400 
700 


Beyond  700  eV  use 
8.75  x  10"15 


-2m  IH2,). 

2 . 2(-16) 

2.5(-16) 

2.94(-16) 

3.5(-16) 

4.9(-16) 

6.04(-16) 

8.2(-16) 

9. 3(-16) 

9.98(-16) 

1. 15(-15) 

1.25(-15) 

1.60(-15) 

1.8(-15) 

1 . 8(— 15) 
1.5(-15) 
1.4(-15) 

1 . 25(— 15) 
1. 10(-15) 
1.05(-15) 
9.8(-16) 
9. 2(-16) 
9.1 ( — 16 ) 
1.0(-15) 
1.05(-15) 
1.0(-15) 
9.4(-16) 
8.4(-16) 
8. 2(-16) 
6.0(-16) 
4. 1(— 16) 
2.0(-16) 
9.0(-17) 
3.5 ( — 17 ) 
1.25(-17) 


